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ABSTRACT 

Spontaneous rapid growth of strong magnetic fields is rather ubiquitous in high-energy density environments 
ranging from astrophysical sources (e.g., gamma-ray bursts and relativistic shocks), to reconnection, to laser- 
plasma interaction laboratory experiments, where they are produced by kinetic streaming instabilities of the 
Weibel type. Relativistic electrons propagating through these sub-Larmor-scale magnetic fields radiate in the 
jitter regime, in which the anisotropy of the magnetic fields and the particle distribution have a strong effect on 
the produced radiation. Here we develop the general theory of jitter radiation, which includes (i) anisotropic 
magnetic fields and electron velocity distributions, (ii) the effects of trapped electrons and (iii) extends the 
description to large deflection angles of radiating particles thus establishing a cross-over between the classical 
jitter and synchrotron regimes. Our results are in remarkable agreement with the radiation spectra obtained 
from particle-in-cell simulations of the classical Weibel instability. Particularly interesting is the onset of the 
field growth, when the transient hard synchrotron- violating spectra are common as a result of the dominant role 
of the trapped population. This effect can serve as a distinct observational signature of the violent field growth 
in astrophysical sources and lab experiments. It is also interesting that a system with small-scale fields tends to 
evolve toward the small-angle jitter regime, which can, under certain conditions, dominate the overall emission 
of a source. 

Subject headings: radiation mechanisms — laboratory astrophysics — high-energy-density physics - ganmia- 
ray bursts — shock waves — reconnection 
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1. INTRODUCTION 

There is a lore that a relativistic particle of charge e and a Lorentz factor 7 moving through a magnetic field B produces 
synchrotron radiation, whose spectrum peaks at cj^ ~ {eBjmc)^^, has an asymptotic ijj^l^ dependence below the peak and falls 
off exponentially at higher frequencies (it makes a second power-law for an isotropic ensemble of particles having a power-law 
distribution in energy). This is often true, but not always. If the field is inhomogeneous on scales comparable or smaller then the 
particle Larmor radius, Ag < 7?^ ~ ^ymc^ /e{B), the produced radiation spectrum may be far different from synchrotron. 

Gamma-ray bursts, supemovae shocks, relativistic pulsar winds and shocks, relativistic jets from quasars and active galactic 
nuclei, magnetic reconnection sites, plasmas produced by high-intensity lasers — they all are the high-energy density envi- 
ronments where conditions are favor able for the spontaneous magnetic field production. The field generati on via the Weibel 



instability ([We ibel 1959; Fried 1959; 



Silva, et aL]|20Q3 Wiersma & Achtenberg 2004 ) or its modifications (Dieckmann et al. 



|2006[ |Bret et al.j20 08 ; Bret 2009 ; Frederiksen & Dieckmann'2008) has been predicted to occur in astrophysical shocks with low 



ambient magnetic field and rare particle collisions, e.g., in gamma-ray burst and large-scale structure shocks ([Medvedev & Loeb 
[1999; Medvedev et al. 2006; Medvedev & Zakutnyaya 2009 ). It has been observed in numerical simulations of relativistic non- 
magnetized shocks ( Nishikawa, et al. 2003^; Frederiksen, et al.^2004 ; Spitkovsk y 2008 ; Keshet et a l. 2009, Nishikawa et al.|20()9] ), 
nonrelativistic shocks unmagnetized and weakly magnetized shocks ( |Kato & T akabe 2008, 2010|, cosmic rays interacting with 
a pre-shock medium fNiemiec et al.^2010|), magnetic rec onnection in electron-positron relativistic and non-relativistic plasmas 
(Zenitani & Hesse 2008 ; Swisdak et al. 2008{ |Liu et al.|2 009), as well as in simulations of and even real laser plasma experiments 
(Ren et al. 2004 ; Tatarakis et al. 2003 ). 

Given such a ubiquity of the process at hand, a natural question to ask is: Are there any observational signatures, which can 
benchmark the process in astrophysical sources and, if any, what can we lear n about the physical conditions there? With the 
radiation t echniques being developed and i mplemented into numerical codes ( |Hededal||2005 [Nishikawa et al.||2Q08| [Sironi & 
Spitkovsky 2009 ; Frederiks en, et al.[[2010| ), we will soon be able to answer this question in detail (for instance, [Medv edev & 



Spitkovsky|200 9 , showed that PIC simulations can realistically model some astrophysical sources). In order to correctly interpret 
the results of simulations and observational data, a comprehensive theory of radiation processes in a strong small-scale magnetic 
turbulence is, therefore, of great demand. 
The e ffects of small-scale inh omogeneities on radi ation emission have been o f long-sta nding theoretical interes t (see. 



for 



instance. 



Landau & Pomeranchuk 1953; Migdal 



1954^ |T956irGinzbu^^^ Syroyatski i 1965; Landau & Lifshitz|197T| and many 
ft Turther been applied to synchrotron radiation from large-scale homogeneous 



more). Techniques developed in these papers have Turther been applied to synchrotron radiation from large- 
magnetic fields wit h a small-scale random field component, as a model of radio emission by cosmic rays in the interstellar 



medium (see, e.g., Nikolaev & Tsytovich| 1979; Bel'kov , et al|1980l [Toptygin'1 985t[Tbptygin & Fleishman|198 7, and references 
therein). The so-called perturbative approach of radiation emission from random small-scale magnetic fields without a large-scale 
component has first been discussed as a model of radiation from Weibel-mediated relativistic colisionless shocks of gamma-ray 
bursts ( |Medvedevj|2000|), where it was referred to as jitter radiation This approach was ge neralized in ( Fleishmann20 06) and 
further correctecP ^ ' ' ' ^ ^ ^ ^ ^ ^ 



([Medvedev"2006 



ie Vi|2000|), where it wa^ 
;^in ^edyedev||2005 
16) was the first to cons 



2006 



Medvedev, 



, et al.||2009||Reynoids et al.||201o) . 
lagnetic turbulence (e.g., Weibel- or 



Because of this and also because 

consider anisotropic magnetic turbulence (e.g., Weibel- or filamentation-instability-generated 
magnetic field tur bulence), we will refer to this paper in the following discussion. We also would like to mention here a recent 
paper by (Revi lle & Kirk|201Q] ), who developed a new algorithm to compute radiation from small-scale turbulent fields. 

In this paper we develop a theory of jitter radiation that accounts for anisotropics of the magnetic field and particle velocity 
distributions, including a trapped population, and further extend the theory to the large angle jitter regime. Our theoretical findings 
are tested with dedicated particle-in-cell simulations. Interesting conclusions are presented in the final section. 

2. THEORY 

Radiation emitted from magnetic fields with small coherence length is not synchrotron, regardless of the actual shape of the 
produced spectrum. In some cases, the spectrum may resemble that of synchrotron, while in others it can be markedly different. 
What kind of spectrum is produced is, in general, set by how curved the particle paths are (i.e., how large their deflections from 
a straight line) compared to the relativistic beaming angle ^1/7. In particular, when the deflecti on angle a ~ e{B)\B/'^mc^ is 



smaller than the beaming angle the particle radiates in the classical jitter regime ( [Medvedev|2000| [2Q06| ), in which the particle's 
velocity v is almost constant, its path is almost straight, r = v^, and its acceleration w = v is random and varies rapidly in time. 
Qualitatively, in the small- angle jitter regime, i.e., when the jitter parameter 



(5j,tt = e{B)\B/mc^ = -/Xb/Rl ^ 5.9 x lO'^ (A^/l cm) {{B)/l gauss) 



(1) 



is small, 6^^ <C 1, the spectrum has a peak at uj ^ (c/Xb)^^- The spectral shape is generally not universal at lower frequencies 
(it can be flat, Fy (xv^, in the isotropic magnetic turbulence but can also be as steep as iv ex i^Mn the presence of a strong 
anisotropy) and is usually a power-law above the peak, whose index is related to that of the spectrum of the magnetic turbulence 
and/or particle distribution. 

2.1. Jitter radiation in the small deflection angle regime 



^ It is important to note here some problems in fFleishman*2006'», which are relevant to the present paper. First, it was argued that the spectrum Ft^ (x. 
below the peak, "valid in the presence of ordered small-scale magnetic field fluctuations, does not occur in the general case of small-scale random magnetic field 
fluctuations." This statement was shown ^Medvedev 1 2005, |2Q0 6 ) to be flawed in that jitter radiation from random magnetic fluctuations with a fairly general 
distribution function (not just "ordered small-scale" fields) does allow for oc spectra. Moreover, the entire range from oc i/^ to oc i/^ is allowed for a single 
electron emission, softer spectra can be expected for some electron energy distributions. Second, another confusing issue is related to the absorption-like oc 1/^ 
jitter spectrum. Such a spectrum is due to plasma dispersion and does not occur in the absence thereof, as the reader might incorrectly infer (see Medvedev 2 005 1 
for more discussion). Finally, we should note that contrary to the claims, none of their papers treat the jitter regime properly. The approximations made in their 
analysis apply to the systems with large-scale plus small-scale isotropic turbulent fields, which do not hold for the generic Weibel-like magnetic turbulence. 
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Here we generalize the theory of classical jitter radiation. Energy emitted by an accelerated relativistic particle and observed 
at infinity is given by t he Poynting flux, which is easily calculated using Lienard-Wiechert (retarded) potentials. One arrives at 
the familiar expression ( [Landau & Lifshitz|1971| ) for the total energy emitted per unit solid angle dO per unit frequency duj: 

dW = ^ l^x[(^-^/^)x 7--]P ^ , (2) 
27rc^ (l-/3v-n)4 2^' 

where I3 = v/c = (1 -7"^)^/^, a "hat" denotes unit vectors, Wc^/ = J we^^ ^ dt is the Fourier component of the particle acceleration, 
the frequency in the comoving and observer's frames are related as cj' = cj(l -/3v • n), and n points toward the observer. Since 
acceleration in a magnetic field, w = {e/^mc)\ x B, is orthogonal to v, the following holds: 

|nx [(n-v/3)xw^.]|' = |w^f (l-/3v-nf-|n-w^f 7-^ (3) 

where we keep the small term ~ 7"^ because it can be important if a particle velocity distribution is structured at angular scales 
< 1/7 and/or if 7 is not very large. 

So far, no approximations were made. Now we use the standard small-deflection approximation, v = const, and we let the 
magnetic field vary, hence the Fourier image of acceleration and its projection onto n are w^, = {ejS /^m)\ea^-^(y^B2j, -v^B^,) 
and ^aW^/ = (el3/'ym)^eaf3'yB'^,(nf3V^-n^Vf3), where eaf3j is the Levi-Civita tensor and in flat configuration space we do not 
distinguish between co- and contra- variant components. Now we immediately obtairj^ 



= {Sa^-VaV^)wf (4) 

|n . Wc^/ 1^ = [Sa(3 (1 - (v • n)^) -VaV(3-nanf3 + (\ • n)(Vanf3+V(3na)] W^P (5) 

where W^P = {ejS /^mfB^fB*^ is the acceleration tensor, B^^/ represents the temporal variation of the field along the particle path 
and is the Kronecker symbol. In general, the field varies in space and time 

Ba,r) = ilTi)-^ j e-'^^'-^ ''^'Rn;^dVid^. (6) 
For a straight path r = ro + v^, the field is a function of one independent variable hence 

B^. = (2^)-4 j e'^'^dt (^-^■^^^-^•^^-^•^o^Bn,k) dQdk 

= (271)-^ J S(uj' -n + k-\)e'^ '''Bn^kdndk, (7) 

where we used the identity for the Dirac (5-function: J e^^^ dt = 27tS(Q. 

So far we dealt with a particular representation of a particle motion through magnetic turbulence. The assumption of ergodicity 
allows us to relate the "representative" history of 5-field along a path to the spatial average over all possible initial positions Fq. 
This is a strong assumption, but it is valid for statistically homogeneous turbulence with no correlation between particles and 
fields (this is not strictly true for trapped particles, see below). Thus 

{K'K^') = (27r)-V-^ J B^^B*^^^^ 8{uj' -Vt + k-y) 8{uj' - l^i +ki • v) e'^^'^'^''' dYQdVtdVtidkdki 

= (27rr^V-^ J B^^Bl^^^S(uj'-n + k-\)dndk, (8) 

where V is the volume occupied by the magnetic field and we again used the identity / e^^^'^^'^'^^dvo = (27r)^(5(k-ki). It is trivial 
to prove via straightforward substitution"^ that the spectral tensor B^ k^rfk is just a Fourier image of the two-point autocorrelation 
tensor of the field, {B^Bi^it^r)) = T'W^ J B^it' y)B(3(t' -t.r' -r)dt'dr' , that is : 

^^,k<k = TvJ (5«5^(^r))^^(^^-^->J^ Jr, (9) 

where T is the duration of an observation. Eqs. Q-pJ and ^ completely determine the jitter radiation spectrum of a single 
particle in the small-deflection regime for an arbitrary neld distribution. 

When magnetic field is generated by a Weibel-type instability, particles bunch into localized streams - current filaments - thus 
determining a preferred direction, s, so that the magnetic field is predominantly in the plane orthogonal to s, as is illustrated in 
Figure [T] so 

^^7,k^?2^k= l^^,k|'(^c./9-W^), (10) 



' Here we used the identities 

^cxX ^afi ^OLh> 

eap-i^x^u = ^/3A (^/3i/ 

^7 A ^7 /J. ^7 



and ecpjecfiu = 



4 Indeed, B^^B^^^ = fBUt'y)e^^^''-^'''^dt'dr'fB*^it''y')e-^^^'''-^^^^ = f[fBo,(t'y)B^(t'-ty-r)dt'dr']e^^^'-^-^^dtdr = 

TV J{BaB^(t,r))e'^^'-^-''^dtdr. 
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s - stream direction 



n - observer direction 



V - particle velocity 



Fig. 1 . — A cartoon showing the assumed planar distribution of the magnetic field typical of the Weibel turbulence and the triad of the unit vectors: s, n, v. 



where \Bq^\^\ is the standard scalar spectrum of the field. The acceleration tensor becomes 



(11) 



where {W^j' \ is the scalar frequency spectrum of the particle acceleration. Finally, the spectral energy emitted by an ensemble of 
particles with a homogeneous anisotropic velocity distribution, F(y), is 



dOduj {lirfc^ 



I 



l+(s-\y 1 (s-vr + (s-nr-2(v-n)(s-v)(s-n) 



(l-/3v-nr T 



7W / (27r 



(l-;9v-nr 
(5(w( 1 - /?v • n) - 1^ + k • v) dVLd\ 



F(\)d\. 



(12) 



The terms in this expression with septuple integration have clear physical meanings. The two terms in the first square brackets 
are the geometric factors coming from the product of (Saf3 -SaS^) in Eq. ( 11 ) with the tensors in front of W^P in Eqs. ^ 

respectively, the term in the second square brackets is simply the particle acceleration spectrum {W^'f' and the outer integral 
weighted with the particle distribution function sums up the contributions of all particles in the system. 

It is worthwhile to note that (i) the radiation spectrum of a single particle is proportional to the spectrum of the particle 
accelerations 

dW/dOdu(x\w^^\^ (13) 

and, hence, is a 'probe' of the magnetic field structure along its path and (ii) the radiation is strongly beamed in the direction of 
the particle motion: 

dW/dO oc (l-/3(n- v))"^ oc (14) 

The latter expression is valid for 7 1 and a small angle ^ between v and the line of sight. If the magnetic field is static (ft = 0) 
and its spatial spectrum has a peak at a characteristic coherence scale, ^ A^^ , then the single electron emissivity in the small 
angle regime is peaked at the frequency 



2.2. Jitter radiation from the trapped population 



(15) 



The Weibel instability ( |WeibeH1959) is driven by anisotropy of the particle distribution F(y), which has been interpreted as 
the instability of streaming particles ( |Fried|1959] ). An infinitesimal transverse modulation of the stream density (i.e., the current 
density) results in transverse magnetic fields, which pinch the particles into filamentary structures and therefore enhance the initial 
perturbation. The filamentary distribution of particle streams is maintained by the self-generated magnetic fields. These stream 
particles are essentially trapped in the filaments for a long time and, hence, the ergodicity assumption used in the derivation of 
jitter radiation in Section [2TT] fails for them, so they are not accounted for by Eq. ( [12| ). These particles oscillate in filaments and 
radiate. Depending on the oscillation amplitude 6m they emit radiation in the small-angle jitter regime if 6^^^ < I/7 and in the 
large-angle jitter regime otherwise. However, in both cases one needs to know exact particle trajectories in order to calculate 
radiation. This is not possible for a generic magnetic field distribution. We illustrate this in the following simplified example of 
a two-dimensional straight filament. 
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If there is translational invariance along the filament, the parallel component of the generalized momentum of a particle is 
conserved: P\\ = 7mv|| -(e/c)A\\(r±) = const, where A\\(r±) is the only non-zero component of the vector potential and it is a 
function of the transverse coordinate only. Since E = in the system, the particle energy is also conserved, therefore v = const. 
Because of the axial symmetry, the magnetic field is zero on the axis r± = 0, so we can take Ay (0) = 0. Then, it is straightforward 
to obtain the equation 

^=[v^-(vii(0)+:^A„(.,)) j , (16) 

which can be solved only in quadrature to obtain the particle trajectory as ^ = t(rj_y, here V||(0) is the parallel velocity on the axis 
of the filament. Even the turning points, r± = r^, can only be found implicitly: v\\(0) + (e/'ym)A\\(rt) = v. 

This model can be simplified further, as shown in Figure [2j For a uniform distribution of current near the filament axis, the 
magnetic field is a linear function of coordinates B(r±) = (r±/ro)Bo, which holds for small displacements rj_ <C tq, where tq is 
the transverse size of a filament and is the field on its "surface". The equation of motion p = (e/c)Y x B for small amplitudes 
and V ~ c reads: 

rl (eBo/jm)(r^/ro). (17) 

This is the equation of a harmonic oscillator. The streaming particles are trapped in filaments and oscillate with the characteristic 
bounce frequency 



fth - {eBo/jmro 



.1/2 



(18) 



Due to the oscillatory motion, the magnetic field in the frame comoving with the particle varies in time as B(t) Bq {vt sin(l]^0 / ro) 
^miVtijt). For small oscillation amplitude^ 6^^^ ~ v^/v <C I/7, the particle path can still be approximated as a straight line, so Eq. 
( 12 ) can be used. To account for the bounce motion, the magnetic field spectrum can be approximately taken in the forrrj^ 




(X 



(19) 



If there is a distribution of filaments of various sizes ro,, and strength Bo^i in the system, the (5-function shall be substituted with the 
actual distribution of bounce frequencies, /^(fi). The observed radiation from the trapped population is peaked at the frequency 



(20) 



2.3. Spectra in the small angle jitter regime 

As we have mentioned earlier, the small-angle jitter spectra are not universal: they depend on the turbulent magnetic field 
spectrum and its anisotropy, in addition to the electron distribution which can also be anisotropic. Figure [3] shows 'Green's func- 
tion' of jitter radiation — the single-electron angle- averaged emission spectrum (equivalent to the ensemble-averaged spectrum 
from monoenergetic randomly moving electrons) from the turbulent magnetic field with a separable ^-function spectrum, i.e., 
|^^7,kP = S(kx-kB)S(ky-kB)S(kz-kB) with kB ~ l/Ag. Other radiation spectra are generally obtained via convolutions with the 
field spectra and electron distributions. For comparison, we also show the single-electron synchrotron spectrum. 

Weibel turbulence is generally anisotropic with the streaming direction and the transverse plane dynamics being decoupled. 
This is confirmed with PIC simulations, see below. The field spectrum can be represented as a separable function |5j7,kP = 
f±(k±)fz(kz) with k\=kl-\-ky. Here we assume f±(k±) ex k~f ^ k± > kB and fz(kz) cx k~^ ^ kz> kB. The spectrum of radiation in 
this case is angle-dependent (see Medvedev 2006; Reyn olds et al. 2010, for a detailed analysis). The spectrum observed along the 
filament direction (z-direction) is shown in Figure [4] with the green sharply -peaked curve resembling the jitter spectrum in Figure 
|3] The low-energy power-law index is 1 and the high-energy index r] is related to the field distribution index sls r] = a. Note that 



Fig. 2. — A cartoon showing a two-dimensional model of a radiating particle trapped in a filament with the magnetic field shown with vectors and the curved 
envelope. The color gradient illustrates the distribution of a current density. The 6r_\_ is the particle oscillation amplitude, the filament size is tq and Bq = B(ro). 



^ This is a more restrictive condition than Srj_ <C ro, for which 6m can also be larger than I/7, depending on the value of tq. 

^ Alternatively, one can set 17 = (static field) and modify the -spectrum of the magnetic field to have a sharp peak at = flh/c. 
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such harder-than- synchrotron uj^ spectrum is also produced by trapped particles regardless of the /^-spectrum. The spectrum of 
radiation observed in the transverse direction (orthogonal to the filaments) is shown with the blue curve. The spectrum is flat 
below the jitter frequency and is a power-law above. The high-energy index is also related to the field spectrum: ^ = /3- 1. The 
spectra observed at intermediate angles are shown with dashed grey curves. 

In the course of non-linear evolution, Weibel turbulence becomes more isotropic. Such 'aged' turbulence is expected, for 
example, in the region far behind a Weibel-mediated shock. The field spectrum becomes isotropic, |5j^,kP = f(k) with = 
kl+ky+k\. We can assume here, for simplicity, that f{k) ex ^ k>kB- In this case, the spectrum is also represented by a broken 
power-law blue curve, which is flat below the break and falls off with the index ^ = /3-2. 



2.4. Jitter radiation in the large deflection angle regime 

When the deflection angle of a particle, 0^ ^ (Sp)/p ^ P±/p^ is greater than I/7 (that is, Sjni > 1), the radiation spectrum is 
determined by the geometry of the particle trajectory. The peak frequency depends on how fast the beaming cone sweeps through 
the line of sight, as in synchrotron radiation. The power of lower frequency harmonics, however, depends on the global structure 
of the path, such as the deflection angle Od or the oscillation amplitude 0^ of a trapped particle. 



The radiation spectrum can be qualitatively obtained as follows (for non-perturbative approach, see, e.g., |Toptyginll985] ). A 
generic particle trajectory can approximately be represented as a smoothly joined set of circular segments (arches) of a certain 
curvature (Larmor) radius, Rl ~ jmc^/e{B), and an angular extent, 6d, as is illustrated in Figure [s] Within each segment, 6 <6d, 
the radiation harmonics are constructed coherently in the way similar to the standard synchrotron radiation. For a harmonics 
constructed over the angle 0, the radiation formation length is / ~ RlO. Radiation emitted over this length is observed in the lab 
frame over the time interval ~ l/cj^, where 7 is the mean Lorentz factor of the particle toward the observer, which is smaller 
than 7 because the transverse motion of a particle is also relativistic: 7^ ~ jO, hence 7^7/7x^1/6^. The observed frequency is 
uje ^ ^ / 0^ c / Rl (jOy^uoj^, where uq = c/Rl is the fundamental (Larmor) frequency. The peak (synchrotron) frequency 
corresponds to the angle 6^^1/7 (i.e., when the observer is within the radiation cone) ujs = uJi/^ ^ (^oj^, hence uje ^ (jOy^uJs- The 



log Fioj) 




0.001 

log OJ 

10"-^ 0.001 0.1 10 

Fig. 3. — Radiation 'Green's functions'. The blue sharp curve is the single-electron angle-averaged emission spectrum (equivalent to the ensemble-averaged 
spectrum from monoenergetic isotropic electrons) from the turbulent magnetic field with a separable 5-function spectrum. The red smooth curve is the single- 
electron radiation spectrum in the synchrotron regime, shown for comparison. The frequencies are normalized to the respective characteristic frequencies, i.e., 
UJj = CJ^ = 1. 
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Fig. 4. — Characteristic small- angle jitter radiation spectra from small-scale magnetic turbulence. The green sharply-peaked curve is the spectrum emitted in 
the direction of the strong anisotropy of the field. This is also the spectrum of the particles trapped in magnetic filaments. The blue double-power-law spectrum 
is emitted in the direction, orthogonal to the anisotropy direction. This is also the spectrum from the isotropic magnetic turbulence. The grey dashed curves 
represent spectra emitted at some intermediate angles with respect to the filaments. 
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radiation spectrum resembles synchrotron spectrum with cj^/^ law below the peak down to the break jitter frequency ujj ~ ^jit^t'^^, 

where Sjut = (jOdT^ by definition. At larger angular scales, > Od, the field is effectively incoherent, therefore radiation occurs 
in the small-angle jitter regime with the effective coherence scale Xb ~ RlOj- 

To summarize this subsection, the jitter radiation in the large angle regime corresponds to intermediate values of the jitter 
parameter: 1 < Sjut < 7. The radiation spectrum is shown in Figure[6] It resembles the synchrotron spectrum, ex uj^^^ exp(-a;/cj^), 
above the break frequency 



{i^b/RlY' {e{B)/mc) 7^ ^ {c/Xb)' {e{B) hmc) 



(21) 



and is similar to the small angle jitter spectrum below this break. At sufficiently large energies, the jitter high-energy power-law 
shows up above the synchrotron exponential decay. Although the small angle jitter spectrum is not universal, it is likely flat ex uj^ 
because the field is generally more isotropic at larger scales and trapping of particles is less efficient, but spectra as hard as cx cj^ 
cannot be excluded and may occur under certain conditions. Note that the break frequency, ujj, depends on the product {B)Xb 
via (5jitt and the spectral peak ujs which measures B (they both proportional to 7^ as well), so one can, in principle, determine the 
magnetic field correlation length Xb from the spectrum alone. Finally, the case of Jjitt > 7 corresponds to synchrotron radiation. 

3. COMPARISON WITH SPECTRA FROM PIC SIMULATIONS 
3.1. PIC setup and results 

The PIC simulations of the relativistic filamentation (Weibel) instability (|We ibel'1959 Fried|1959| have been performed and 
radiation from the system has simultaneously been obtained (see, |Frederiksen, e t al. 2010| for technical details of the simulations). 
The simulations represent the classical Weibel instability with two equal density charge-neutral electron-positron and electron- 
ion plasma streams. Both two-dimensional and three dimensional setups with various initial bulk Lorentz factors, F, in the range 
of 2 to 15 were use d. In three-dimens ional electron-positron pair plasma simulations reported here (an exhaustive description is 
presented by |Frederiksen, et al. 201Q| the simulation box of 500^ cells with periodic boundary conditions and with the resolution 




Fig. 5. — Cartoon showing the piecewise circular decomposition of a particle trajectory. For each segment of a certain curvature (Larmor) radius, Rl, and the 
deflection angle, 6^, radiation is produced at times between and + ^^- Radiation is beamed directly toward an observer on the right within a highlighted 
~ 1/7 part of the segment. 




O.Ol 



log OJ 

Fig. 6. — Characteristic radiation spectra in the large-angle jitter regime. The spectrum is synchrotron-like above the jitter break and resembles the small-angle 
jitter spectra below the break, cf. Figure|4] Note that the jitter high-energy power-law can show up above the synchrotron exponential decay at sufficiently large 
energies. The jitter break frequency is uoj = S'^^^^s, where (5jitt > 1 is the jitter parameter and ous is the synchrotron frequency. 
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Fig. 7. — The magnetic field distribution in 2D simulations, for illustration purpose only. The simulation box is 2000 x 40 skin lengths, snapshots are at times 
r ~ 2, 4, 6, 8, 10, 12 ujpl. The field magnitude (with the polarity: into and out of the plane) is coded with white and yellow colors. The pair of numbers at the 
bottom of each panel shows min and max values of the field in simulation units. Saturation occurs around f ~ 12 in this particular run; before this the filaments 
grow in amplitude but not in transverse size. Note that at saturation, the local field amplitude keeps increasing for some time due to the nonlinear evolution and 
pinching of the filaments. The overall field energy density is decreasing, however, due to the decreasing filling factor. 

of 10 cells per the relativistic skin length, 5e = ujpe/c = {Ane^n/TmeC^) , was used. All physical spatial and temporal scales are 
expressed in units of the skin length and the plasma time, cj~] , respectively. Spectra are collected in situ during runtime using the 
PhotonPlasma code developed at the Niels Bohr Institute (Haugb0lle 2005 ; | Hededal|2Q05| . In the simulations reported here, 
the total number of particles in the simulation box was ~ 10^^; however, the spectra were collected from about N ~ 10^ particles 
and sub-cycling A^rad < 10~^ A^pic with linear trajectory interpolation was used to resolve high radiation frequencies. The spectra 
were calculated coherently from the retarded electric fields of the ensemble of emitting particles 



duo df] 



/ ^E,,,je^^^'dt^ 



(22) 



where the phase is (j)' = t' - ii ■ ro(t')/c and we neglected self-absorption, Rasin and other plasma effects. Thus, the obtained 
spectra represent 'time-resolved' radiation emitted within the time interval {to^ti}. Snapshots of the system (particles, fields, 
radiation) are saved every two plasma times. Here we report the 'head-on' case, when the radiation is emitted along the initial 
streaming direction. The results for oblique and edge-on cases are reported elsewhere (Frederiksen, et al. 2010). 

The PIC simulations cover the initial exponential growth of the magnetic field, saturation of the Weibel insta bility and the 
nonlinear evolution and mergers of current filaments (the merger model is developed in Medvedev, et al. 2005 ; Shvets et al.|20Q9 



when the magnetic field gradually decays. The distribution of the magnetic filaments during the instability growth and saturation 
is illustrated in Figure |7] (this is the illustration only; the conditions of the scientific runs are different, see above). The evolution 
of the field strength at aifferent stages of the instability is shown in Figure [8] where 

eB = B^/ [87rr(r - l)nmeC^] , (23) 

the magnetic field energy normalized by the total initial kinetic energy in the system, is plotted. The field strength reaches 
maximum ^ 10% at about saturation, ^ ~ 10. After saturation, the field strength averaged over the simulation box decreases as 
B(t) (X . Note, however, that the local field amplitude can still increase for some time after saturation due to the nonlinear 
evolution and pinching of the filaments, but the overall field energy density falls because of the decreasing filling factor. 

The particle distribution function, i^(7Vx,7V||), and the transverse and parallel (with respect to the streaming direction) mag- 
netic field spectra \Bk^ p = {B_\_(k±)B*^(k±)) and l^^^u p = {B _\_(k\\)B*^(k\\)) are shown in Figure^^at representative times: during 
the linear stage of the instability with exponential growth (^ = 6), at saturation (t = 10) and thelate nonlinear regime (t = 30). 
Hereafter, we consider F = 10 case only; other cases with relativistic F are very similar. 

The particle distribution is well described by a multi- Gaussian in jVxj jVy^ jv^. It is represented by two shifted Gaussians in 
the streaming (parallel) direction and a Gaussian in the perpendicular plane. The parallel and transverse temperatures are different 
and vary in time, so does the average (bulk) Lorentz factor of the interpenetrating streams. 
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Fig. 8. — The magnetic field energy density normalized by the initial kinetic energy density of the streams, cb, versus time in realistic 3D simulations for two 
runs with T = 2 (blue) and F = 10 (red). At late times, oc r^. 

Concerning the electromagnetic field, the electric fields and parallel magnetic fields are present but very weak (typically, less 
than a percent), so we only use in calculations of the field spectra. It was also found that the field spectrum is separable, that 
is, it can be expressed as a product of two independent functions: |5kP = \Bk^ Pl^fcy ^^^^ being evaluated independently. In 
order to account for the entire evolution, we used a model of a Gaussian plus a smoothly broken power-law. Since the number of 
the fit parameters is rather large, some of them (but different) were kept fixed at various times. At early times, the perpendicular 
spectra, \Bk^ p, are fit well by a Gaussian or a double-power-law with a broad transition region whose width quickly reduces 
with time. At saturation, it is still a double-power-law. After saturation, the peak of the spectrum is much better described by a 
Gaussian component, which moves toward lower-/:, leaving behind a single power-law at higher k. At late times, the power-law 
exhibits some curvature (downturn) at small k, where it intersects with the Gaussian component. In the beginning, the parallel 
spectra, \Bk^^ p, are flat (nearly at the noise level) with a broad peak at small k and then they develop a power-law. This peak is 
broad and is always at the largest scale (smallest k) corresponding to the box size. The peak is caused by the "causality horizon": 
the filaments at distances greater than the light crossing time are uncorrected thus producing white noise at small k. The peak 
width is decreasing as ~ 1/^, therefore. The bump at large k is likely due to the numerical Cherenkov instability. The power of 
these Fourier modes is very small and they do not significantly affect radiation spectra (although some signatures of it can be 
discerned). Overall, except for the linear phase of the instability, the temporal evolution of the parallel spectrum is very modest. 

3.2. PIC and semi-analytic radiation spectra 



Semi-analysical small-angle jitter spectra are obtained from Eq. ([12]) for each PIC snapshot using the available data for the 
particles F(\) = F(jv±,jv\\) and fields l^k,^^!^ = f±(k±)f\\(k\\)fb(^), where f±(k±) = \Bk^\^, /||(^||) = l^^n P and the bounce 
frequency distribution is taken heuristically, using Eq. ([18]), as follows. If all the filaments are identical, then the total transverse 
size of the filament ~ 2ro corresponds to the peak of the f±(kj_) distribution, 2ro ~ k']^. We conjecture that the profile of fj_, 
especially near the peak, describes the ensemble of filaments in the system. This is certainly violated at kj_ greater than the skin 
scale beyond which no filaments exist. These scales are suppressed with an exponential factor in our model, therefore. Taking 

into account that Q - c/ro - 2c/(2ro) - 2k^c, we have Mft) = [2nf^(2n)e-^^/'^pT(r-l)/j^] ^'"^ . Although not rigorous, this 
prescription works well, likely because the main effect is due to the presence of a characteristic frequency whereas the exact 
shape of fh^S^) is of lesser importance. We have also computed the synchrotron spectra for the same plasma parameters, i.e., the 
spectra one would expect from the large-scale magnetic field of the same strength with the same electron energy distribution. 
These spectra are similar to the large-angle jitter spectra, except for the low-frequency part, where a shallower spectrum is 
expected in the jitter regime (see discussion in Section [24]) 



The comparison of the PIC and semi-analytic spectra are shown in Fig. 10 for two time intervals: t G {4, 12} which corresponds 
to the early exponential growth before saturation and r G {14,26} which corresponds to the late nonlinear phase of the filament 
merger and field decay. The predicted small-angle jitter and classical synchrotron spectra are plotted for comparison. The early- 
time PIC spectrum is in agreement with the semi-analytic prediction but not with the synchrotron spectrum both in the position 
of the peak and the overall spectral shape. In particular, (i) the PIC spectrum exhibits synchrotron-violating ~ J^-^ scaling law 
below the peak with the subsequent flattening at lower frequencies, (ii) the high-frequency part in the PIC spectrum does not 
show the near-exponential roll-off (note that the electron distribution is still close to monoenergetic at this time), and (iii) the 
width of the peak region in the PIC is substantially narrower than it is in the synchrotron spectrum. The late-time PIC spectrum 
is generally consistent with synchrotron, except at frequencies below cj ~ 10, where some flattening becomes evident 

The fact that the late-time PIC spectrum is synchrotron-like isn't surprising because (5jitt > 1 in the run (see Section 2.4 ), which 
is evident from Figure [TT] where we plotted the jitter parameter 



^jitt 



;[2e5r(r-l)]i/2(/:skmA±,max), 



(24) 



where the peak of the -spectrum, ^±,max, is found numerically at each time and ^skin is the wave-number corresponding to 
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Fig. 9. — A representative sample of snapshots of (top to bottom) the particle distribution, /^(7v^,7V||), the perpendicular, \B]^^ p, and parallel, \B]^^^ 

magnetic field spectra before, at and after saturation (left to right: t = 6, 10, 30) for the 3D PIC simulation with T = 10. The blue curves are the empirical models, 
while the red dots are computed directly from the data. 

the skin scale. What is particularly interesting is that the early-time PIC spectrum is not consistent with synchrotron at all. The 
explanation to this is simple. The particle distribution is highly anisotropic at these times: there are still well resolved particle 
streams which form and are also trapped in the current filaments. For these particles, the size of the magnetic structures is 
irrelevant. As long as their angular deflections are small enough, they emit in the small-angle regime and produce hard spectra, 
no matter how large the ^jitt parameter is. We have confirmed it by using tracer particles in a "frozen field" simulation. The full 
snapshot of particles and fields at r = 8 was used and in the subsequent simulation, the particles were allowed to move, but the 
fields were kept fixed (i.e., static and, of course, not self-consistent anymore). The evolution was traced over fourteen plasma 
times - long enough for the particles to become substantially isotropized, so that the number of trapped particles diminished 
drastically. The radiation spectrum obtained in this run is entirely consistent with the synchrotron spectrum from the particles 
with a thermal spread, as is evident from Figure [12] 

4. CONCLUSIONS 

The primary results of this paper are as follows. First, we present a general expression for the spectral energy per solid angle 
emitted by an ensemble of particles in the small-angle jitter approximation, Eq. ([12]). Second, we h ave fou nd that the electrons 
streaming through the filaments and being trapped in them produce a transient ham spectrum, Eqs. ( |18|20| ). Third, we analyzed 
the large-angle deflection regime and showed that the spectrum starts to resemble the synchrotron spectrum near the peak, but a 
new spectral break at a lower frequency appears, Eq. pT) . From the positions of the spectral peak and break, one can deduce the 
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Fig. 10. — Radiation spectra at early, t E {4, 12}, (left panel) and late, t G {14,26}, times (right panel). The spectra obtained 'on-the-flight' from 3D PIC 
r = 10 simulations of the Weibel instability in the e^e'-pair plasma are the two 'noisy' blue curves. The semi-analytic spectra from Eq.[T2]are shown in red and 
the synchrotron spectra are in green. The early time PIC spectra are consistent with the hard (synchrotron violating) small-angle jitter emission and the late time 
emission is consistent with the large-angle jitter or synchrotron. 




10^ FT 




Fig. 12. — The PIC spectrum from the "frozen field" simulation (dark blue) and the original spectrum from Fig.Qofleft) (light blue) are shown (in both, F = 10). 
The magnetic field structure is nearly identical in both, but the trapped population is practically absent in the "frozen-field" case. The corresponding synchrotron 
spectrum for the actual particle distribution (green) and for the monoenergetic electrons (dashed grey) are also shown for comparison. 



field correlation length. 

Fourth, PIC simulations show that the radiation spectrum produced at the onset of and during the phase of the exponential 
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growth of the magnetic field is grossly inconsistent with synchrotron, Fig. \T0\ The appearance of such a spectrum in the 
beginning of an emission episode can be used as a benchmark signal of the onsetof the magnetic field generation in astrophysical 
sources and laboratory experiments. Among possible astrophysical systems where such emission can be or could have already 
been obser ved are gamma-ray bursts. The data show hard synchrotron- violating spectra in some bursts and the majority of spectra 
are flat (Preec e et al.|1998||2000[[Kaneko et al.'20Q6), which are difficult to explain within the synchrotron model. Such spectra 
have recently been interpreted in the jitter emission paradigm ( Medvedev, et al. 2009 ). We should caution the reader that "blind" 
application of the jitter spectrum template for interpretation of observational data, without checking the physical conditions at 
and the validity of the jitter approximation for the source in hand (galactic and quasar jets, supernovae remnants, etc.), can yield 
incorrect results. Possible laboratory experiments include laser-pl asma interactions in which a beam (e.g., a probe electron beam) 
propagating through turbulent fields can emit jitter radiation ( Reynol ds & Medvedev|201l] ). 

Fifth, although the spectra after saturation are consistent with the synchrotron 'template , an overall trend of the system toward 
the small-angle jitter regime (i.e., toward Jjitt < 1) is observed in Fig. UU which suggest the scaling: Sjut oc ^/e^r(ujpt)~^-^ . 
Although the magnetic field decays, it does so rather slowly, B(t) ex (ujpt)\see Fig. [8](for example, for F = 10 the systen shall 
return to the small-angle jitter regime at times t > 100). So, if the field is continuously produced (as in the case of the propagating 
shock, for instance) the field decay can be compensated by the increase of the emitting volume, so that the total spectral emissivity 
will increase logarithmically, Ptoti^) oc J BdV oc ln(ujpt) oc ln(ujpLps/c), assuming these scalings hold at asymptotically late times, 
where Lps is the size (longitutinal extent) of the post-shock medium. Since the plasma time ^ uj~^ is generally very short in 
astrophysical sources, the overall time-integrated (and even time resolved, but with a coarse temporal resolution) spectrum can 
be dominated by the small-angle jitter spectrum, which is expected to be flat (unless the strong anisotropy of the fields and/or 
particles is somehow maintained). 

Here we also comment on the relation of our results and the simulations of (ISironi & Spitkov sky 12009'). First of all, one should 
understand that the physical set-ups are entirely different. In this paper we study the evolution of radiation during the filamentation 
instability: the plasma in the simulatio n box is initially h omogeneous, the particle distribution is unstable and the entire system is 
intrinsically non- stationary. In contrast, |Sironi & Spitkovsky ( 2009 ) simulated a well-developed collisionless shock: a steady state 
system with very slow, if any, evolution of conditions and radiation spectra. Thus, the two simulations are complementary and 
the direct comparison of them is not well-posed. Since, however, the shock is moving through a medium with a constant speed, 
the temporal evolution of the filamentation instability and its subsequent saturation and further nonlinear evolution studied here 
is, to a certain degree, represented by the spatial profiles in the pre-shock and post-shock domains. The filamentation instability 
and its saturation occur far in front of the shock (hundreds or thousands skin lengths, in typical simulations) by particles escaping 
from the shock, hence this region would roughly correspond to the early and saturation times in our simulations. After saturation, 
mergers of magnetic filaments increase their sizes and radiation spectrum now mimics synchrotron, as our simulations show. In 
shock simulations, this merging stage occurs in a large region in front of the shock and this is where the radiation is collected in 
( fSironi & Spitkovsky ,2009j ). Their results are, thus, in agreement with the ones presented here. The shock itself and the medium 
just behind the shock do not correspond to our simulations. However, the magnetic field strength decreases behind the shock and 
in a few hundred skin lengths radiation should enter the jitter regime. Unfortunately, radiation from neither the early pre-shock, 
nor from the far downstream regions have been shown in (Sironi & Spitkovsky 2009). We stress that radiation from the far 
downstream can be of great importance and dominate the entire shock emission, provided the scalings presented in the previous 
paragraph hold. 
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